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INTRODUCTION 

When  a  plane  shock  wave  moving  through  one  material  meets  an  interface 
with  a  new  material,  the  shock  is  transmitted  at  a  different  strength  into 
the  "target"  material  and  a  new  wave  is  reflected  back  into  the  "parent" 
material.  The  strength  of  the  transmitted  wave  is  determined  in  this  report 
for  the  case  when  both  materials  are  solid  and  can  be  described  by  a  Mie- 
Gruneisen  equation  of  scate. 

The  reflected  wave  may  be  either  a  shock  or  a  rarefaction.  In  the 
first  case,  the  transmitted  shock  pressure  exceeds  that  in  the  incident 
shock.  This  occurs  generally,  but  not  necessarily,  when  the  target  material 
is  denser  than  the  parent  material.  In  the  second  case,  the  reflected 
rarefaction,  there  is  a  reduction  in  shock  pressure.  This  situation  arises 
typically,  but  not  necessarily,  when  the  target  material  is  of  lower  density 
than  the  parent  material.  One  result  of  this  report  is  to  exhibit  the 
boundary  between  these  two  cases.  This  boundary  can  also  be  interpreted 
as  the  relation  between  bulk  modulus  and  density  in  the  target  material 
such  that  there  is  no  reflected  wave  and  consequently  all  of  the  energy  in 
the  incident  wave  is  transmitted. 

This  analysis  was  made  to  facilitate  estimates  of  the  damage  incurred 
in  structures  which  are  strongly  shocked,  as  from  hypervelocity  projectile 
impact,  intensive  x-ray  exposure,  or  underground  nuclear  explosions.  Actual 
failure  of  structural  materials  is  incurred  as  a  result  of  shock-wave 
reflection  at  either  a  free  surface  or  at  a  contact  surface.  The  resulting 
negative  pressure  could  lead  to  tensile  fracture  and  possibly  spallation 
fragments.  However,  failure  criteria  will  not  be  considered  in  this  report. 
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REFLECTION  OF  ACOUSTIC  WAVES 

When  a  wave  is  sufficiently  weak,  whether  tensile  or  compressive, 
it  is  referred  to  as  an  acoustic  wave  and  its  analysis  can  be  carried  out 
in  considerable  detail  since  the  governing  equations  are  linear.  The 
general  solution  to  the  acoustic  reflection  problem  is  discussed  separately 
in  this  section  since  an  elementary  analytical  solution  is  available. 

The  solution  to  the  nonlinear  shock  reflection  problem  discussed  in  the 
body  of  this  report  has  the  acoustic  solution  of  Eq.(l)as  a  limit  when 
the  shock  pressures  are  well  below  the  bulk  modulus  of  the  material.  For 
the  geometry  outlined  in  Fig.  1,  the  ratio  of  the  pressure  after  reflection, 
p  ,  to  the  pressure  before  reflection,  p^,  is 
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where  p  is  the  density,  c  is  the  sound  speed  based  on  the  bulk  modulus, 
c  =  a/k/p  ,  K  is  the  bulk  modulus,  and  pc  =  ^Kp  is  the  acoustic  impedance. 
The  result  given  in  Eq.  (1)  shows  that  the  pressure  ratio  never  exceeds 
2  in  the  acoustic  case. 

The  above  result  can  be  used  to  determine  the  pressure  ratio  for  a 
wave  transmitted  through  a  slab  of  material  embedded  in  a  different 
material.  If  the  subscripts  1  and  2  are  used  to  denote  the  outside  and 
inside  materials,  respectively,  then  the  ratio  of  transmitted  pressure, 
P^,  to  incident  pressure,  P^,  is 
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NOTE:  A  zero  subscript,  (  )  ,  denotes  the  state  of  unshocked  parent 
material;  although  no  such  material  appears  in  the  figure,  it 
will  have  bulk  modulus  KQ  and  density  p  . 


Fig.  l--Geometry  for  the  shock  reflection  problem. 


If  the  acoustic  impedance  ratio  is  as  large  as  10,  or  as  small  as  1/10,  the 
pressure  transmitted  through  the  sandwich  is  64$  of  the  original.  This  is 
about  the  greatest  reduction  i/hat  might  be  expected  with  common  materials. 
The  effect  of  subsequent  reflections  within  the  sandwich  may  complicate  the 
result  for  some  applications. 
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REFLECTION  OF  SHOCKS  IN  IDEAL  GASES 

For  ideal  gases,  the  px’essure  ratio  behind  the  reflected  shock  for 
the  geometry  of  Fig,  1  is  given  by  Courant  and  Friedrichs ^  where  the 
target  material  is  treated  as  a  rigid  wall; 
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where  p^  is  the  initial  static  pressure  and  |i  =(y-l)/(y+l)  depends  on 
y,  the  ratio  of  specific  heats.  For  strong  shocks,  this  simplifies  to 


p3  "  po  3y  -  l 
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00 


which  is  8  for  y  =  7/5  and  6  for  y  =  5/3.  For  weak  shocks,  Pq/P^  i8  near 
unity  and  the  pressure  ratio  is  near  2,  as  in  the  acoustic  case  when  the 
target  material  has  infinite  acoustic  impedance.  Thus,  the  shock-ref leC'/.on 
pressure  ratios  are  larger  in  the  ideal  gas — strong  shock  case  than  in  the 
acoustic  case.  It  is  this  fact  that  motivated  the  work  described  here, 
since  for  strong  shocks  in  solids  it  was  not  clear  whether  the  large  pressure 
ratios  found  in  the  perfect  gas  case  or  the  small  pressure  ratios  found  in 
the  acoustic  approximation  would  be  more  nearly  correct. 
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EQUATIONS  FOR  THE  REFLECTED  SHOCK  CASE 

The  geometry  for  the  general  shock  reflection  problem  is  outlined  in 
Fig.  1,  which  defines  the  four  regions  that  appear  in  the  shock-reflection 
prob.em,  (Material  in  the  initial  state,  denoted  by  a  zero  subscript,  is 
not  sh'wn  in  Fig.  1.)  The  usual  shock  relations  as  given  in  Ref.  1,  page  129, 
are: 


Pi  v2  =  p2  v2 


P1+„1V1  =P2+P2v2  , 


Pi/Pi  +  \  +  I  \  =  P2/P2  +  \  +  %  v2  ^ 

where  denotes  the  relative  velocity  of  the  material  and  the  shock.  In 
addition  to  these  conservation  equations,  an  equation  of  state  of  the  Mie- 
Gruneisen  form, 

p  =  G(V)  E/V  +  F(V)  ,  (8) 

is  assumed  in  order  to  complete  the  physical  description  of  the  problem. 

A  review  of  the  experimental  data  in  Ref.  2  suggests  that  if  the  Grunelsen 
ratio  G(V)  has  the  form 

G(V)  =  Gq  V/VQ  ,  (9) 

then  the  data  are  in  many  instances  fitted  better  than  if  G(V)  is  taken  to 
be  a  constant.  This  assumption  also  leads  to  convenient  analytical  ex¬ 
pressions  in  the  calculation  of  isentropes.  In  the  present  analysis,  GQ 
is  taken  to  be  2,  which  is  typical,  but  not  necessarily  accurate,  for  most 
materials.  The  best  choice  of  GQ  for  various  materials  generally  lies  in 
the  range  from  1.5  to  3.0,  but  the  final  results  appear  insensitive  to  the 
value  of  G^  as  indicated  in  Fig.  2. 

(2) 

The  determination  of  F(V)  follows  from  the  empirical  observation' 
that  the  shock  velocity,  ug,  and  the  material  velocity,  u^,  behind  a  shock 
into  a  material  at  rest  are  related  through  the  equation 


u  =  e  +  8  u  . 
s  p 


(10) 
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pj_gt  2 _ The  effect  of  variation  in  the  constant  G~  on  shock  reflection  pressures 

for  equal,  hulk  modulus,  Kp  =  K_,  and  three  values  of  density  ratio. 
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The  empirical  parameter,  S,  usually  has  values  in  the  range  from  1.2  to  2.0 
and  was  taken  to  he  1.5  in  the  calculations .  As  with  G^,  subsequent  results 
appear  insensitive  to  the  value  of  S,  as  illustrated  in  Fig.  3.  Combining 
this  equation  with  the  shock  relations  of  Eqs.  (5),  (6),  and  (j)  leads  to 

f(v)  =  kb(i  -  gq  e/2)/(i  -  S9 )“  =  k  f(e)  ,  e  =  i  -  v/v  (11) 


where  K  is  the  bulk  modulus  of  the  material  and  V  is  the  initial  specific 
volume.  Equations  for  the  pressure  and  velocjty  in  Region  4  behind  the 
transmitted  shock  can  be  obtained  from  the  first  two  shock  equations,  (5) 
and  (6),  in  terms  of  the  compression  0^,  where 

04  =  1  -  VV2  •  <12) 

(The  initial  specific  volume,  V,  is  in  this  instance.)  These  relations 
are; 


f(V 

P4  =  *2  1  -  Gq  6^/2 

and 

\  s  /0  4  P4  V2  * 

(it  will  be  assumed  that  the  Gruneisen  ratio  is  the  same  in  both  materials.) 

A  comparison  of  the  shock  pressure  from  Eq.  (13)  with  the  more  detailed 
treatment  by  Tillotson  in  Ref.  3  is  given  in  Fig.  4,  which  shows  the 
Eugoniot  for  pressures  up  to  100  Mbar  in  copper.  (Copper  is  a  convenient 
example,  since  it  has  the  se2.ecte<3  values  of  G^  =  2  and  S  =  1.5.) 

Equations  for  the  compression  0^  and  the  pressure  p^  in  Region  3  are 
obtained  by  considering  the  conservation  equations  across  the  reflected  shock. 
The  appropriate  equations  relating  these  conditions  in  Regions  3  and  1  are, 
after  appropriate  manipulations  of  Eqs.  (5)  through  (9), 
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Fig.  4 — Comparison  of  shock  pressure  in  the  low  pressure  approximation  of 
Eqs,  11  and  12  with  the  more  general  Hugoniot  of  Ref.  3. 
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where 


i  -  o0(u3  -  ex)/2 

e3  -  1  -  Vvo  '  8i  -  1  '  Vvo 


(16) 


(17) 


The  additional  requirement  that  the  pressure  and  velocity  be 
continuous  between  Regions  3  and  4  completes  the  formulation  of  the 
problem. 
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EQUATIONS  FOR  THE  REFLECTED  RAREFACTION  CASE 

Region  3  of  Fig.  1  may  contain  a  rarefaction  wave  under  some 
conditions.  If  this  is  the  case,  the  pressure  will  vary  isentroplcally 
in  that  region,  and  it  is  necessary  to  have  an  equation  relating  pressure 
and  specific  volume  on  the  isentrope.  This  is  obtained  by  considering 
the  equation  of  state 

P  =  G0  E/Vq  +  F(V)  .  (18) 

In  order  that  dS  =  (dE  +  pdV)/T  be  a  perfect  differential,  it  is  necessary 
that 


3  A  9E\  _  Ji_ 
av  'T  3T'  "  3T 


1  ,  3ES 

t  (p  + 


or,  substituting  from  the  equation  of  state, 

-  +  T  GqE^Vq  =  G  E/V0  +  F(V)  , 

where  the  subscripts  denote  differentiation  with  respect  to  that  variable. 

Ih) 

The  characteristic  equations'  '  are 


.  dV  _  y  dT 


dE 


rr  =  v0  TG0  =  F(V)  +  G0  E/V0  • 

Two  independent  solutions  containing  arbitrary  constants  c7  and  c0  are 

x  cL 

obtained  by  combining  the  first  equation  with  the  second  and  then  that 
equation  with  the  third; 


GqV/V, 

c  =  E(V)  +  E  e  u  ,  where  H(V) 


,  fV  ,( 


F(V)  dV  , 


°2 =  °o  n +  ^  t  * 

leading  to.  the  general  solution 
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where  J  denotes  an  arbitraiy  function.  H(V)  is  plotted  in  nondimensional 
form  in  ^ig.  5  for  the  function  given  by  Eq.  (11).  If  the  specific  heat 
is  assumed  constant  and  equal  to  Cy,  then 


5E 

ST 


=  j'(T 


Gov/vo 


) 


and  therefore  the  function  J  must  be  linear  in  its  argument.  It  follows 
that  the  equation  of  state  is 


E  =  Cy  T  ~  e 


and  the  entropy  is  given  by 


°oV/Vo  /  e°oV/V°  F(v) 


dV 


(20) 
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(21) 


where  SQ  is  an  arbitrary  constant.  Combining  this  expression  with  the 
equation  of  state,  the  isentrope  connecting  the  pressures  in  Regions  1  and  ^ 
is  given  by 

[Pl  -  F(V1)]eG°Vl/V°  +  Gq  H  (V1)/VQ  =  [p3  -  F^)^0^  0 


+  G0  H<V/V0 


(22) 


(1) 


The  Riemann  integral'  '  is  used  to  find  the  relation  between  u  and  V 
along  the  particle  paths.  The  appropriate  solution  for  the  current  problem 
is 


ux  -  u3  -.J  SJla  W(V1)  -  W(V3)  . 


(23) 


The  sound  speed,  c,  appearing  in  the  above  integrand  can  be  obtained  from 
the  equation 
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(24) 


and  the  expression  below  for  the  pressure  on  an  isentrope  from  the  point 
(V,  p^)  in  the  p-V  plane; 

-G  V/Vn  ,  G  V,/V  , 

p  =  F(V)  +  e  0  0  {gq  [H(V)  -  H(V3)]/V0  +  lp3  -  F^)]  e  0  3  °}  . 


(25) 
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CALCULATIONS 

The  equations  for  shock  transmission  were  solved  hy  an  iteration 
technique  using  the  IBM-7040.  A  general  outline  of  the  solution  procedure 
is  given  in  the  flow  diagrams  on  the  two  following  pages.  The  calculations 
were  done  using  nondimens ional  pressures  and  velocities,  which  are  denoted  by 
bars.  The  following  list  relates  the  quantities  in  the  flow  diagram  to  those 
used  in  the  theoretical  discussions  above. 

i>i  =  p-j/Kq  p2  =  p2/k2  p  =  p3/kq  ,  =  p4/k2 

\  =  V°0  ^  =  Vc2  S  =  U3/C0  '  *4  =  VC2 

c0  =  K(/P0  C2  "  K2^P2 
1(0)  =  H(V)/K0  ,  6  =  1  -  V/VQ  . 

W(0)  =  ~-W(V) 

C0 


SOLUTION 


Chart  1.  Solution  procedure  for  the  reflected  shock 
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Chart  2,  Solution  procedure  for  the  reflected  rarefaction. 
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SUMMARY  AND  CONCLUSIONS 


The  purpose  of  this  report  is  to  provide  a  survey  of  the  pressure  ratios, 
p^/pp  encountered  in  the  shock  reflection  process.  In  Figs.  6-12  the  results 
of  the  calculations  are  given  for  initial  shock  pressures  in  the  range  from 
zero  to  ten  times  the  bulk  modulus,  with  both  the  density  and  bulk  modulus 
ratios  being  varied  in  the  two  materials.  When  the  transmitted  pressure  is 
less  than  that  in  the  incident  wave,  the  reflected  wave  is  a  rarefaction,  and 
when  it  exceeds  that  in  the  incident  wave,  the  reflected  wave  is  a  shock.  The 
curves  for  the  two  cases  fit  smoothly  together,  as  might  be  expected  on 
theoretical  grounds,  by  noting  tliat  when  the  pressure  ratio  is  near  unity  the 
reflected  wave  is  weak,  and  thus  the  Hugoniot  and  isentropes  for  the  reflected 
wave  have  a  third-order  contact. ^  The  pressure  ratios  lie  between  1/3  and 
3  in  the  region  of  the  parameters  Investigated,  which  includes  bulk  modulus  and 
density  ratios  from  1/10  to  10.  The  velocity  ratios,  u^/u^,  are  plotted  in 
Figs.  13  through  19. 

The  "complete  transmission"  curves  separating  the  reflected  shock  case 
from  the  reflected  rarefaction  case  are  shown  in  Fig.  20  as  a  relation  between 
bulk  modulus  ratio,  K2/Kq,  and  density  ratio,  Pg/P0*  This  relation  depends 
on  shock  strength.  When  the  shock  is  weak,  the  complete  transmission  curve 
is  given  by  acoustic  theory  and  is  KgP2  =  KQp0.  As  the  incident  shock  strength 
is  increased,  the  curve  separating  the  two  cases  approach  the  line  p2/pQ  = 

Thus  for  strong  shocks  complete  transmission  requires  only  that  the  density 
ratio  be  near  unity.  The  density  ratio  is  therefore  more  important  than  the 
bulk  modulus  in  determining  the  transmitted  shock  pressures  for  strong  shocks. 
This  is  to  be  expected  on  physical  grounds,  since  at  higher  pressures  materials 
begin  to  behave  as  a  gas  and  the  bulk-modulus  terms  in  the  Mie-Gruneisen 
equation  become  relatively  small. 


Fig.  9--The  ratio  of  transmitted  pressure,  p^,  to  incident  pressure 
the  case  where  the  hulk  modulus  ratio,  Kg /%->  equals  1. 


transmitted,  pressure,  p_,  to  incident  pre&sure 
■e  the  hulk  modulus  ratio,  K^/Kq  equals  2. 
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Fig.  18 — The  ratio  of  particle  velocity  in  the  transmitted  wave,  u_,  to  particle  velocity  in 
the  incident  wave,  u_ ,  in  the  case  where  the  bulk  modulus  ratio,  Ko/Kh,  equals  h. 
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Fig.  19 — The  ratio  of  particle  velocity  in  the  transmitted  wave,  u_,  to  particle  velocity  in 
the  incident  wave,  u^,  in  the  case  where  the  bulk  modulus  ratio,  K^/K^  equals  10. 
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Fig.  20 — The  boundary  between  the  reflected  shock  case  and  the  reflected 
rarefaction  case  as  a  function  of  bulk  modulus  ratio,  Kp/K-.,  and  density 
ratio,  Po/Pq*  5116  ko’-odazy  depends  on  shock  strength,  p^,  which  is  ex¬ 
pressed  in  units  of  the  bulk  modulus,  K0,  in  the  diagram.  If  the  para¬ 
meters  are  such  that  a  point  lies  on  the  boundary  there  is  no  reflected 
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